Enhanced squeezing with parity kicks 
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Using exponential quadratic operators, we present a general framework for studying the exact 
dynamics of system-bath interaction in which the Hamiltonian is described by the quadratic form of 
bosonic operators. To demonstrate the versatility of the approach, we study how the environment 
affects the squeezing of quadrature components of the system. We further propose that the squeezing 
can be enhanced when parity kicks are applied to the system. 



Introduction - Coupling between system and environ- 
ment is ubiquitous in all quantum processes (e.g., in 
quantum information processing). Such coupling usually 
results in: (i) the energy decay of the quantum system; 
(ii) the destruction of the relative phases of several super- 
posed quantum states, and thus the linear superposition 
of several quantum states turn into a classical mixture. 
However, the environment can also help us, for example, 
entanglement between two systems can be generated via 
a common environment [l[ . 

Although it seems impossible to model the environ- 
ment exactly in many cases, and thus difficult to ob- 
tain the exact dynamics of a system-environment inter- 
action, quantitative analysis based on approximate de- 
scription of the environment is needed in many cases, for 
instance, the analysis of decoherence suppressing meth- 
ods H 0, i, i, @, 0, H IE Gil ; the discussion of the entan- 
glement of two systems coupled to the environment; and 
the study of the quantum dissipation of systems |ll|, ll2| . 
An extensively adopted approach to model the environ- 
ment, which is also called a reservoir or bath, is to in- 
troduce a set of harmonic oscillators with different fre- 
quencies. In this case, the interaction between the sys- 
tem and the environment is modeled by coupling the sys- 
tem to these harmonic oscillators through an appropri- 
ate interaction Hamiltonian. Several methods have been 
proposed to study the coupling between the system and 
a set of harmonic oscillators. In quantum optics (e.g., 
Refs. [HI, [3, GI||), a quite often used method to ana- 
lyze Markovian process is either a master equation or 
a Langevin equation. Another method is the path in- 
tegral approach (lrij which was extensively developed 



in Refs. [ll|, [17[, but, this method is very complicated. 



Furthermore, different approximations are used in all of 
these methods to make the problem tractable for either 
analytical or numerical calculations. 

In this paper, we introduce a new method to calcu- 
late the evolution of the bosonic system, coupled to the 
environment. The total Hamiltonian is described by a 



quadratic form of the bosonic operators. Our method is 
based on some properties of exponential quadratic oper- 
ators. As shown in below, this new method provides a 
feasible way to calculate the effect of the environment on 
system. As an example, we apply our method to study 
the environment effect on the generation of squeezed 
states. Moreover, we also use our method to study the 
system-environment interaction when the parity kicks are 
applied to the system. We find that the parity kicks can 
help us to obtain a better squeezing. 

Exponential quadratic operators. — For a set of anni- 
hilation operator a^(l < i <jn), exponential quadratic 
operators(EQO) (see 18, ljj 2Q|) are expressions of the 
form 
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Here [ai,aj] = 1. The above equation can also be 



e *A ra j n w J 1 j c J 1 



written in the following way Q 
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ric matrix. If we define S 



then we have 



QA T Q- 1 = A 



T e -RS 



(2) 



where the multiplication in eQA T e~® is understood to 
act on each term of A. 

Coupling between oscillator and reservoir. — Consider 
a system comprising of a harmonic oscillator with anni- 
hilation operator a, and a reservoir consisting of a set of 
oscillators with annihilation operator b k for each mode. 
The Hamiltonian of system-reservoir is described by 

H = huata + ^2 fakblh + h^ lk {ab\ + b k o)), (3) 

k 
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where the first, second, and third terms are the system, 
reservoir, and system-reservoir interaction Hamiltonians, 
respectively. Here, ^ k are the coefficients representing the 
coupling strength between the system and the mode k of 
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reservoir - these coupling constants are typically much 
smaller than the other frequencies in the Hamiltonian. 
For simplicity, but without loss of generality, we regard 
these couplings as reals. 

We calculate the evolution of a{a*) in the Heisen- 
berg picture by using equation @. For U = e ~ lHt l h , 
putting A T = (a^, b\ , b\, ■ ■ ■ , 6^, 0,61,62, ■ • ■ , b n ) and 
P^ 

P 



frame). 



R = 



where 



P 
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(4) 



Thus, to calculate the evolution of a (a'), we need only 
to calculate the matrix e~ RS . 

Coupling between system and reservoir during a 
squeezing process. — To see the power of the technique, 
let us consider a Hamiltonian for degenerate parametric 
amplification with a classical pump under the influence 
of a reservoir in a squeezing process. The Hamiltonian 
can be expressed as 



H = Wa + Uhe[e 2iu)t a 2 - e^Vf 1 



fkOjb^bj + h{o) 7j6j + h.c). 
3=1 j=i 



(5) 



In order to remove the time dependence in the Hamilto- 
nian, we transfer the Hamiltonian into a rotating refer- 
ence frame with U — cxp(iHot/h) with Hq = h\o(a) a + 

n 

b\bj). Thus in the rotating reference frame, the 

3=1 

Hamiltonian in Eq. ((5]) becomes 

1 ™ 
Hi = - -ihe[(a^) 2 - a 2 } + ^ h(u % - u)b]bj 



3 = 1 



(6) 



K^^ljbj + h.c). 



Note that the first term is just the usual squeezing op- 
erator. We can easily find the matrix R corresponding 
to —iHit/H. By analyzing e~ RS , numerically if necessar- 
ily, we obtain the evolution of a*{t) and a(t), and thus 
the solution of all quantities associated with a squeez- 
ing process. The most important one among them is 
((A(aW+a(«)t)) 2 } = ((AI) 2 }. 

Parity kicks in the squeezing process. — Using appro- 
priate time varying control fields, it is well known that 
one could alleviate decoherence effects through a se- 
quence of frequent parity kicks. As in Ref. [J], we in- 
troduce an extra Hamiltonian (in the rotating reference 



Hi — Hi + Hkick (t) , 



where H kick (t) = H klck for U < t < U + t and 
H k i ck = otherwise. We require tj+i — t, = To for all 
i. Moreover, we assume r <C To and i?feicfe(i) is strong 
enough during the kick periods that we can neglect the 



effect of Hi, which is e 



-iHir/h 



iH kiek r/H Undcr 



these conditions, we will model parity kicks as unitary 
operators P = e - lH kickr/h ac ^ m g on S y S tem at a set 
of time ti. Since we want to eliminate the influence 
of coupling between system and reservoir, we require 
P to have following properties PH system P = H system , 
PH bath P = H bath , and PH int P = -H int . The three 
Hamiltonians are defined in Eq. (J5]). It is easy 
to verify that P = e - l7laia satisfies above equa- 
tions. Thus the unitary operator corresponding to 
two such periods would be p e -' lH '^/ h p e -^ H '^/ h = 



D -(iTo/h)(H svstc 



i+H bat h,-Hint) e -(.iT /h)(H S y Bte 



Y. Intuitively, it shows the interaction between system 
and reservoir of different periods cancel each other out. 
In fact, it has been proved that when tq — > the system 
and the reservoir are totally decoupled. 

We use numerical computation to verify this effect in 
the squeezing process. To this end we calculate the evo- 
lution of (v(t) and a(t): We have 

o t (2nT ) =y +n o t y n . 

To use the EQO method shown in Eq. ([2|) to solve the 
above expression, we note that if 



e yi A T e-^ 



A T Pi, 



then 



e Y2 e yi A T e- yi e- r2 =A T P 2 P 1 . 

Thus, we know that we need only to calculate the e~ RS 
in 

Y^A T Y = A T e~ RS , 

and e~ nRS would be the desired transforming matrix. 
Again we use the above property and see that we only 
need to calculate 

e (ir /h)(H systcm + H bath ±H int ) A T e -(iT /h)(H systcm +H bath ±H int ) 

(7) ' 

For simplicity we consider the ground state situation. 
The procedure is entirely general and applies for the 
case of T > OK. We compute the variance {(AX) 2 ) 
with two types of coupling: namely, the Lorcntzian spec- 
trum and the ohmic spectrum. For the Lorcntzian spec- 
trum 7j = g(u)j) = rjT / (ujj — u>) 2 + T 2 , as an example 
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FIG. 1: The variance ((AX) 2 ) versus rescaled time et when 
parity kicks are on and off respectively. Figure (a) refers to a 
Lorentzian spectrum and figure (b) refers to an Ohmic spec- 
trum. We find that a better squeezing can be obtained when 
parity kicks are applied to the system. 



of numerical calculations, we assume T = 2 x 10 9 Hz, 
i] = 5 x 10 7 Hz, the squeezing parameter e = 10 8 Hz and 
the kick period tq = 1.67 x 10~ 9 s. For the ohmic spec- 
trum 7j = g(ujj) = v / ^cUje~^/" c , we assume £ = 10 6 Hz, 
u> c = 10 9 Hz, the squeezing parameter e = 7x 10 8 Hz and 
the kick period tq = 2.5 x 10 _9 s. For both spectrum, 
we assume the frequencies associated with the system 
and reservoir to be u = 10 9 Hz and uij = j X 10 7 Hz 
(j = 1,2,... 200) respectively. With these parameters, 
the variance ((AX) 2 ) versus rescaled time et is plotted 
in Fig. [TJ which shows that a better squeezing can be 
obtained if parity kicks are applied to the system. 

Discussions. — It is interesting to do a comparison be- 
tween several methods, including the widely used Marko- 
vian master equation (l3| . Under the Hamiltonian Q 



with a Lorentzian spectrum as shown in Ref. 2l| , w 
obtain an exact solution 

a(t) =[u{t)e- rt / 2 a + %(*)^]e" iw * 

={[cos(9;/2) + | sm(6i/2)] + £ u^b^e 



(8) 



The constant is given by 9 = ^An^D - T 2 , where D 
is the density of reservoir modes and Uj (t) are some com- 
plicated functions 21 1 . However, the master equation, in 
the Markovian approximation, can be written as 

^ = - iuj[a)a, p] + ^ [2apa f - a) ap - pa ! a}+ 
+ \n[a) pa + apa) — a*ap — paa)], 



FIG. 2: (a) The probability of the system being in state |1) 
with a Lorentzian spectrum. We can see that our numerical 
solution is close to the exact solution. On the other hand, we 
can see master equation is not valid for this situation. (We 
set the reservoir to have 200 equally distributed oscillators 
while getting the numerical solution of our method.) (b) The 
probability of the system being in state |1) when the coupling 
strengths ■jk are constant. We can see that the lines of our 
numerical solution and master equation's solution coincide 
with each other. 



where A = 2nDg(u>) 2 is a constant, which represents 
the decay rate of the harmonic oscillator. For simplic- 
ity, we assume the temperature of reservoir to be zero, 
and the initial state of system to be |1). We then cal- 
culate the the probability P(t) = (l\trn(p(t))\l), which 
can be used to observe the decay of system. We can 
also compute P(t) when the coupling strengths 7fc are 
constants. For example, we assume that parameters of 
the Lorentzian spectrum in Fig^a) to be 7j = g(ujj) = 
2.8209 x lO^/V^- - u) 2 + 10 12 Hz, and the flat spec- 
trum in Fig H;b) to be 7, = 5.6419 x 10 6 Hz with 
j = 1,2,... 200. For the flat spectrum we assume the 
frequencies associated with the system and reservoir to 
be uj = 10 9 Hz and ujj = j x 10 7 Hz (j = 1,2,... 200) 
respectively For the Lorentzian spectrum, however, 
we change the frequencies of the reservoir to be u>j = 
(50 + j/2) x 10 7 Hz (j = 1,2,... 200) due to the shape 
of the Lorentzian spectrum, which varies dramatically at 
the center and is negligible at two sides. By making this 
change we can sample the spectrum better. Then we 
plot Fig. [2j We can find that while the master equation 
leads to a good approximate solution in some cases, it 
fails sometimes. Thus our method is more reliable, and 
the accuracy can be further improved by using better 
numerical methods. 
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We also note that the parity kicks can be done by in- 
creasing the frequency of the harmonic oscillator for a 
short time interval. For example, this can be achieved 
in the ion trap by changing the electric field, (also see 
Refs. (22], [23| for schemes of generating squeezed states 
in ion trap) 

Conclusions. — We have shown that for a general 
Hamiltonian with bononic quadratic forms, we can com- 
pute dynamics of system using exponential quadratic op- 
erators. Our method provides substantial improvement 
over computation involving master equations as we do 
not need to solve any differential equations and provides 
numerical solution for Hamiltonians that can be written 
in quadratic form of creation and annihilation operators. 
Thus, this new techniques compares well with the dy- 
namics of the system under a master equation but it is 
in some sense more appealing as it could provide in prin- 
ciple analytical expressions for some cases. In particular, 
we analyze the effect of reservoir in a squeezing process 
and we propose possible scheme to improve the degree of 
squeezing. Our method can be applied to study the prob- 
lem on the quantization of nano-mechanical systems, the 
further work will be presented elsewhere. 
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